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Abstract 

In this paper, the impact of energy constraints on a two-hop network with a source, a relay and a 
destination under random medium access is studied. A collision channel with erasures is considered, and 
the source and the relay nodes have energy harvesting capabilities and an unlimited battery to store the 
harvested energy. Additionally, the source and the relay node have external traffic arrivals and the relay 
forwards a fraction of the source node’s traffic to the destination; the cooperation is performed at the 
network level. An inner and an outer bound of the stability region for a given transmission probability 
vector are obtained. Then, the closure of the inner and the outer bound is obtained separately and they 
turn out to be identical. This work is not only a step in connecting information theory and networking, 
by studying the maximum stable throughput region metric but also it taps the relatively unexplored and 
important domain of energy harvesting and assesses the effect of that on this important measure. 
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I. Introduction 


Taking advantage of renewable energy resourees from the environment, also known as energy 
harvesting, enables unattended operability of infrastmeture-less wireless networks. There are 
various forms of energy that ean be harvested, ineluding thermal, solar, aeoustie, wind, and even 
ambient radio power Q. Energy harvesting is reeently seen as a promising feature for wireless 
networks regarding self-sustainability and also effieieney. This permits long-term operation of 
distributed wireless eommunieation systems, sueh as sensor networks, without the need for regu¬ 
lar maintenanee. However, the additional funetionality of energy harvesting in wireless networks 
introduees several ehanges and ealls for assessment of the system long-term performanee sueh 
as in terms of the throughput and stability. The ideal seenario is to make the energy limitations 
transparent to the network. 

Among distributed eommunieation protoeols, we are partieularly interested in ALOHA, a 
simple random aeeess seheme in whieh transmission attempts are performed randomly, indepen¬ 
dently, and distributively Q. In 0, the eapability of energy harvesting was first introdueed in 
the analysis of the slotted ALOHA for a simple setting as an initial step to understand its impaet 
on the aehievable stability region. Reeently, this result has been generalized in [[^ by taking into 
aeeount the multi-paeket reeeption eapability at the reeeiver and finite eapaeity batteries at the 
energy harvesting sourees. In Q, a eognitive aeeess protoeol was studied for the seenario where 
the higher priority primary souree is powered by harvesting energy whereas the lower priority 
seeondary souree is assumed to have a reliable power supply. 

Cooperative eommunieation is one of key teehnologies to aehieve eoverage extension and 
throughput enhaneement in wireless networks Q. In this work, we eonsider paeket-level eoop- 
eration rather at the physieal layer, in whieh a relay node takes responsibility of paeket delivery 
for those it eould overhear and sueeessfully deeode from the transmissions by souree node |(^- 


[10|. A key differenee between physieal-layer and network-layer eooperation is that the latter 
ean eapture the bursty nature of traffie. The impaet of network-level eooperation in an energy 
harvesting network with a pure relay (without its own traffie) under seheduled aeeess (time 
division multiple aeeess in a eontrolled manner) was studied in pT|. 
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A major limitation of Information Theory is the inability to handle bursty traffie and queuing 
delay. In the communications networks bursty traffic and delay are central and indispensable 
concepts. The information theoretic capacity region is derived under the assumption of saturated 
queues. However, under stochastic and bursty traffic arrivals, the maximum stable throughput 
or stability region becomes a meaningful and relevant measure of rates in packets per slot in 
wireless networks. Thus, the maximum stable throughput is an important performance measure, 
akin to information theoretic capacity, but simpler to track and analyze and more appropriate 
for systems with sources that generate signals randomly in time. Understanding the relationship 
between information-theoretic capacity and stability region has received considerable attention 
in recent years and some progress has been made primarily for multiple access channels p^ . 

The characterization of random access stability for bursty traffic is a challenging problem 


even without energy harvesting [13|-[15|. Additionally, for a network with more than three 
users (interacting queues), the exact characterization of the stability region is not known. This 
is because each node transmits and, thereby, interferes with the others only when its queue 
is non-empty. Such queues are said to be interacting with each other in the sense that the 
service process of one depends on the status of the others. The analysis for the case with energy 
harvesting becomes significantly more challenging because the service process of a node depends 
not only on the status of its own queue and battery, but also on the status of the other node’s 
queue and battery. 

In this paper, we study the impact of energy constraints on a two-hop network with a source, a 
relay and a destination under random medium access as shown in Fig. We assume a collision 
channel with erasures. Both the source and the relay node have external traffic arrivals. The 
relay forwards a fraction of the source node’s traffic to the destination and the cooperation is 
performed at the network level. In addition, both source and relay nodes have energy harvesting 
capabilities and an unlimited battery to store the harvested energy. We provide necessary and 
sufficient conditions for the stability of the considered network as shown in Fig. We first obtain 
an inner and an outer bound of the stability region for a given transmission probability vector. We 
then take the closure of the inner and the outer bound separately over all feasible transmission 
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probability vectors. Interestingly, it turns out that the bounds are tight in terms of the closure, 
as also stated in 0 This study provides insights on designing a relay-assisted network under 
energy constraints. When the aggregate charging rate is above one and the source and the relay 
lie in the intermediate traffic regime, the system has identical performance with that of a network 
without energy constraints, meaning that in that regime the energy limitations are transparent to 
the network operation. In this paper we focus on a simple network, as mentioned earlier, more 
realistic and complex systems are impossible to analyze, primarily due to the difficulty in tracking 
interacting queue. However, insights can still be obtained, even from simple models. This work 
provides a step in connecting information theory and networking, by studying the maximum 
stable throughput region metric. Further, it taps the relatively unexplored and important domain 
of energy harvesting and assesses the effect of that on this important measure. 

The rest of this paper is organized as follows. In Section we define the stability region, 
describe the channel model, and explain the packet arrival and energy harvesting models. In 


Section IIT we present inner and outer bounds on the stability region as well as the closure of 
the stability region. The proofs of our results are given in and |Vj Finally, we conclude our 
work in Section I 


II. System Model 

We consider a time-slotted system in which the nodes randomly access a common receiver 
and both source and relay nodes are powered from randomly time-varying renewable energy 
sources, as shown in Fig. Each node stores the harvested energy in a battery of unlimited 
capacity. We denote with S, R, and D, the source, the relay and the destination, respectively. 
Packet traffic originates from both S and R, and because of the wireless broadcast nature, R 
may receive some of the packets transmitted from S, which in turn can be relayed to D. The 
packets from S that fail to be received by D but are successfully received by R are relayed by 
R. A half-duplex constraint is imposed here, i.e. R can overhear S only when it is idle. 

Each node has an infinite size buffer for storing incoming packets and the transmission of 
each packet occupies one time slot. Node R has separate queues for the exogenous arrivals 
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and the endogenous arrivals being relayed through R. Nevertheless, we ean let R have a single 
queue and merge all arrivals into a single queue as the achievable stable throughput region is 


not affected [16|. This is due to the fact that the link quality between R and D is independent 
of which packet is selected for transmission. 

The packet arrival and energy harvesting processes at S and R are assumed to be Bernoulli 
with rates A 5 , 5s and Xr, 5r, respectively, and are independent of each other. Qi and Bi, i = S,R, 
denote the steady state number of packets and energy units in the queue and the energy source at 
node i, respectively. Furthermore, a node i is called active if both its packet queue and its battery 
are nonempty at the same time, which is denoted by the event At = {{Bi 7 ^ 0 } fi {Qi ^ 0 }} 
and idle otherwise (denoted by Ai). In each time slot, nodes S and R attempt to transmit with 
probabilities qs and qr, respectively, whenever they are active. Decisions on transmission are 
made independently among the nodes and each transmission consumes one energy unit. We 
assume a collision channel with erasures in which if both S and R transmit at the same time 
slot, a collision occurs and both transmissions fail. The probability that a packet transmitted by 
node i is successfully decoded at node j{^ i) is denoted by which is the probability that 
the signal-to-noise ratio (SNR) over the specified link exceeds a certain threshold for successful 
decoding. These erasure/outage probabilities capture the effect of random fading at the physical 
layer. The probabilities psD, Prd, and psR denote the success probabilities over the link S — D, 
R — D, and S — R, respectively. We also assume that node R has a better channel to D than S, 


i.e. Prd > Psd- 

The cooperation is performed at the protocol (network) level as follows: when S transmits a 
packet, if D decodes it successfully, it sends an ACK and the packet exits the network; if D 
fails to decode the packet but R does, then R sends an ACK and takes over the responsibility 
of delivering the packet to D by placing it in its queue. If neither D nor R decode (or if R 
does not store the packet), the packet remains in S^s queue for retransmission. The ACKs are 
assumed to be error-free, instantaneous, and broadcasted to all relevant nodes. 
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Fig. 1: A relay-aided wireless network with energy harvesting eapabilities. 


The average serviee rate for the souree node is given by 


/is = [gs(l - gR)Pr {Bs 7^ 0, Ar) + qsPT{Bs 7 ^ 0, x \psD + (1 - Psd)psr] , (1) 

and for the relay is given by 


Pr — [qr^^ — 5 's)Pr {Br 7 ^ 0 , ^ 5 ) + qRPr^BR 7 ^ 0 , ^ 5 )] x prr. ( 2 ) 


Denote by Q\ the length of queue i at the beginning of time slot t. Based on the definition 


in [151, the queue is said to be stable if 


lim Pr[Q\ < x] = F{x) and lim F(x) = 1 

t^OO X^OO 

Loynes’ theorem [ [T7| states that if the arrival and serviee proeesses of a queue are strietly 
jointly stationary and the average arrival rate is less than the average serviee rate, then the queue 
is stable. If the average arrival rate is greater than the average serviee rate, then the queue is 
unstable and the value of Q\ approaehes infinity almost surely. The stability region of the system 
is defined as the set of arrival rate veetors A = (Ai, A 2 ) for whieh the queues in the system are 
stable. 
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III. Main Results 

This section presents the stability conditions of a network consisting of a source and a relay 
both having energy harvesting capabilities, and a destination, as depicted in Fig. The source 
and the relay are assumed to have infinite size queues to store the harvested energy. 

The next proposition presents an inner bound on the stability region by providing sufficient 
conditions for stability. 

Proposition III.l. If (A5, Xr) G TZinner, with 


TZinner = {(As, Xr) : As < min (^s, qs) [1 - min {5r, qr)] [pso + (1 - Psd)psr] , 

^ ^ (1 - Psd)psr 

PSD + (1 - Psd)psr 

then the network in Fig. is stable. 


As < min {5r, qr) [1 - min ((5s, qs)] Prd f , (3) 


Proof: The proof is given in Section IV-A ■ 

The following proposition describes an outer bound of the stability region by obtaining 
necessary conditions for stability. 

Proposition III.2. If the network in Fig. A stable then (As, Xr) G IZ, where TZ = TZi 1J7^2, 
with 


TZi — I (As, Xr) : 
min((5s, qs) \psD + (1 - Psd)psr] 


1 - t - 


min((5s, gs)(l - Psd)psr 


As+ 


[1 - min (( 5 s , gs )] Pi?D 
(1 - Psd)psr 


[1 - mm{6s,qs)]pRD 
Xr < min((5s, qs) [psd + (1 - Psd)psr] , 


Xr + 


PSD + (1 - Psd)psr 


As < min(5jj, qr ) [1 - min(55, gs)] Prd \ ■ 


(4) 









As < min((5R, qR)pRD, 


7^2 


I (As, Ah) : Ah + 


[1 - min((5H, Qr)] (1 - Psd)psr + min((5H, qR)pRD 
[1 - min((5H, gn)] [psd + (1 - Psd)psr] 


As < min(5s, qs) [1 - min(5H, gn)] \psd + (1 - Psd)psr]} 

(5) 


Proof: The proof is given in Seetion IV-B 


Fig. |2(a)| and |2(b)| illustrate the TZi and 7^2 deseribed in Proposition |IIL2 


In the previous proposition we provided the stability eonditions for given transmission prob¬ 
abilities qs and gn, this stability region is denoted by £(gs, gn, Note that 


Thinner ^{qS,qR,^S,^R) ( 6 ) 

The elosure of the stability region is defined by 


^{Ss,Sr)= C{qs,qR,5s,SR). (7) 

(gs>te)e[o,i]^ 

The following theorem deseribes the elosure of the stability region for the network we eonsider. 

Theorem III.l. If Ss + Sr > 1, the closure of the stability region, £(Ss, Sr), is illustrated in 
Fig. 0 and is described by three parts, (i) The line segment AB, where xa = 0, i/a = SrPrr, 
and xb = {1- Sr)'^ [psn + (1 - Psd)psr], Vb = S%prd - (1 - - Psd)psr- 

(ii) the curve from B to C which is described by 


As 


+ 


PsrO- ~ Psd)ISs 


+ 


\R 


P, 


= 1 


( 8 ) 


RD 


Psd + (1 — Psd)Psr y Prd [Psd + (1 — Psd)Psr] 

(Hi) the line segment CD where xc = Sg [psD + (1 - Psd)psr\, 2/c = (1 - Ss)‘^Prd - <^1(1 “ 

PSDJPSR ana xr mm {i-psd)psr ’ {i-Ss)prd+Ss{i-psd)psr /’ 

If Ss + Sr < 1, the closure of the stability region, £(5s, Sr), is illustrated in Fig. and 

is described by the line segments EF and FG, where xr = 0, i/e = SrPrd, xp = (^s(l — 
Sr) [psd + (1 — Psd)psr]> Uf = ^r(1 — Ss)prd — <5s(l — ^r)(1 ~ PsD)psRy 


xg = mm 


il-Ss)SRPRD[psD + {'^-PSD)PSR] ^sH-^s)PRd\PSD + {^-PSd)PSr] 1 _ n 

— PSd'JPSR ^ n n-kAo n —T 1 c n ^ Tl o it> I i/Cj 


il-Ss}PRD+Ss(^-PSD)PSR j 
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(a) TZi 



Fig. 2: An outer bound of the stability region IZ 


'TZi\J'Tl 2 , deseribed in Proposition 


IIL2 


Proof: The proof is given in Seetion |Vj ■ 

Remark 1: Regarding the stability region £(^ 5 , we only have inner and outer 


bounds and not the exaet expression, however for the elosure C{5s,Sr), we do have the exaet 
eharaeterization. 
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Fig. 3: The closure of the stability region for 63 + Sr > 1. 


Remark 2: If ^5 + < 1 as depicted in Fig. the closure of the stability region has a linear 

behavior, which is an indication that the performance of the system is affected by the low energy 
harvesting rates. Furthermore when ^5 + 5^^ > 1, when the arrival rate at the source or the relay 
is at high arrival rate regime then the performance is affected by the energy harvesting rate and 
is depicted by the linear segments in Fig. The interesting case is when both the arrival rates 
Xs and Xr lie in the intermediate arrival rate regime, then the performance is identical to the 
relay network without energy limitations and the closure of the stability region has a non-linear 
behavior. 


IV. Analysis 

To derive the stability condition for the queue in the relay node, we need to calculate the total 
arrival rate. There are two independent arrival processes at the relay: the exogenous traffic with 
arrival rate Xr and the endogenous traffic from S. Denote by Sa the event that S transmits a 
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packet and the paeket leaves the queue, then 

Pr(S'A) = [1 - [psD + (1 - Psd)psr] ■ (9) 

Among the paekets that depart from the queue of S, some will exit the network because they 
are deeoded by the destination direetly, and some will be relayed by R. Denote hy Sb the event 
that the transmitted paeket from S will be relayed from R, then 

Pt{Sb) = [1 - gRPriAR)] (1 - Psd)psr- (10) 


The eonditional probability that a transmitted paeket from S is relayed by R given that the 
transmitted paeket exits node S^s queue is given by 


Pt{Sb\Sa) = 


(1 -Psd)psr 


PSD + {^-Psd)psr 

The arrival rate from the source to the relay is 


( 11 ) 


^S^R — Pi{Sb\Sa)\s. 


( 12 ) 
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Note that the average arrival rate from the source to the relay, Xs^r, does not depend on 
the transmission probabilities or the energy harvesting rates, it depends only from the success 
probabilities psD and psR. 

The total arrival rate at the relay node is given by 

, I i^-PSD)pSR , 

^R,total — ^R-\ - —p. -^ 

PSD + (1 - PSDjPSR 

A. Sufficient Conditions 

A queue is considered saturated if in each time slot there is always a packet to transmit, i.e. 
the queue is never empty. Assuming saturated queues for the source and the relay node, the 
saturated throughput for the source node is given by 

Ps = ~ {Bs 7 ^ 0, Br 7 ^ 0) + gsPr {Bg 7 ^ 0, B^ = 0)} [psn + (1 ~ Psd)psr\ ■, 

(14) 

and for the relay is given by 

Pr = -qsW^{.Bs ^0 ,Br^ 0) + g/jPr(5s = 0, 5^? 7 ^ 0)}prd. (15) 


Each node transmits with probability qi, i = S, R, whenever its battery is not empty and each 
transmission demands one energy packet. Each energy queue i is then decoupled and forms a 
discrete-time M/M/1 queue with input rate 5i and service rate qi, thus the probability the energy 
queue to be empty is given by 


Pr(5,^0) 



(16) 


Note that the battery queues are decoupled since we assume that each node, either the source 
or the relay, attempts to transmit independently in a random access manner. 

Then, after some calculations, we obtain that the saturated throughput for the source is 


/i| = min (^ 5 , qs) [1 - min {5r, qn)] [psD + (1 - Psd)Psr] , 


(17) 
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and for the relay is 


/i^ = min ((5 r, Qr) [1 - min (( 55 , qs)] Prd- 


(18) 


The suffieient eonditions (TZinner) for the stability are obtained by A 5 < /i| and XR^totai < Pr 
and are given by Q, in Proposition |IIL1[ 

The saturated throughput that we obtained in this subseetion is an inner bound of the stability 
region, sinee the assumption that the souree and the relay have always paekets to transmit leads 
to lower aehievable rates in terms of paeket per slot [ [T^ , p9| . 


B. Necessary Conditions 

The average serviee rates for the souree and the relay are given by Q and Q, respeetively. 
The average serviee rate of eaeh queue depends on the status of its own energy and also the queue 
size and the energy statues of the other queues. This eoupling between the queues (both paeket 
and energy) results in a four dimensional Markov ehain whieh makes the analysis eumbersome. 


Therefore, the stoehastie dominant teehnique [1^ is essential in order to deeouple the interaetion 
between the queues, and thus to eharaeterize the stability region. Thus, we first eonstruet parallel 
dominant systems in whieh one of the nodes transmits dummy paekets when its paeket queue 
is empty. Note that even in the dominant system a node eannot transmit if the energy souree is 
empty (beeause even the dummy paeket eonsumes one energy unit). 

We eonsider the first hypothetieal system in whieh the souree node transmits dummy paekets 
when its queue is empty and all the other assumptions remain intaet. The average serviee rate 
for the relay given by Q beeomes 


Pr = {?i?(l -qs)'^^{Bs y^0,BRy^ 0) + qRPT^Bs = 0 ,Br^ 0)}prd. 


(19) 


The average serviee rate of the relay, pR, in the first hypothetieal system is the same with the 


saturated throughput of the relay obtained in (18). From Loyne’s eriterion, the relay is stable if 

Xr, total < PR, thus 


Xr + 


(1 - Psd)psr 
PSD + {1-Psd)psr 


Xs < min {5r, qR) [1 - min ((^ 5 , qs)] Prd- 


( 20 ) 
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The average number of packets per active slot for i? is [1 — min { 63 , ^ 5 )] QrPrd, thus the 
fraction of active slots is given by 

{^-Psd)psr 


Pr {Bh ^ ^ 0) = 


Ad -I- y'^-PSDiPSR 
^ PSD + {i-PSD)PSR ^ 


[1-min {6s,qs)]qRPRD' 

After changing ( [T^ and (21) into ([^, the service rate for the source becomes 

A -I- i^-PSD)PSR \ ' 

^ PSD + ii-PSD)PSR ^ 


ps = min ( 55 , qs) 


1 - 


[PsD + (1 - Psd)psr] ■ 


[1-min {6s,qs)]PRD 

The queue in S is stable if A 5 < ps and after some manipulations we obtain 


1 + 


min {6s, qs) (1 - Psd)psr 
[1- min {6s, qs)]PRD 


R 


^ min ((5g, qs) \psD + (1 - Psd)psr\ ^ 

^ [1-min {6 s, qs)] Prd 

< min (^ 5 , qs) [psD + (1 - Psd)Psr] ■ 


( 21 ) 


( 22 ) 


(23) 


The derived stability conditions from the first hypothetical system are summarized in (|^. 

In the second hypothetical system, the relay node transmits dummy packets and all the other 
assumptions remain intact. Thus, the average service rate for the source given by ([^ becomes 


Ps = {<?s(l - gi?)Pr {Bs ^0,Biij^0) + qsPT{Bs j^0,Bji = 0)} [pso + (1 - Psd)psr] , 


which is equal to saturated throughput of the source and is given by 


(24) 


Ps = min (^ 5 , qs) [1 - min {6it, qn)] [psn + (1 - Psd)psr] • (^5) 

From Loyne’s theorem, the queue in source is stable if \s < ps, thus 

As < min ((5s, qs) [1 - min {6r, g^)] [psd + (1 - Psd)psr\ ■ (^ 6 ) 

The average number of packets per active slot for S is qs [1 — min {6r, g^)] [psd + (1 — Psd)psr]- 
The fraction of active slots for the source S is 


Pr(i?s^0,gs^0) 


_As_ 

qs [1 - min {6r, gs)] \psD + (1 - Psd)psr] ' 


(27) 
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After replacing from ( fT^ and into Q, the service rate for the relay is 

As 


liR = min {6r, qR) 


[1 - min {5r, qr)] [psn + (1 - Psd)psr] 


Prd- 


(28) 


The queue in the relay node R is stable if \r, total < Pr and after some manipulations we 
obtain 

[l-min((5R,gij)](l-psD)psi? + min((5R,gij)pi?D, ^ \ 

Ar +-^-TTT-^^-As < mm [bR, qr) prd- (29) 

[1 - mm (()/?, qR)\ \psD + (1 - Psd)Psr\ 

The derived stability conditions from the second hypothetical system are given by Q. 


An important observation made in [ 141 is that the stability conditions obtained by using the 
stochastic dominance technique are not merely sufficient conditions for the stability of the original 
system but are sufficient and necessary conditions. However, the indistinguishability argument 
does not apply to our problem. In a system with batteries, the dummy packet transmissions affect 
the dynamics of the batteries. For example, there are instants when a node is no more able to 
transmit in the hypothetical system because of the lack of energy, while it is able to transmit in 
the original system, thus it may result to a better chance of success for the other node. 

The obtained stability conditions are necessary conditions of the original system and are 


summarized in Proposition III.2 


V. Proof of Theorem HIT II 

In this section we will derive the closure of the outer of the stability region defined in 
Proposition III.2[ The closure will be obtained over all the feasible transmission probability 
vectors (g^, g^) G [0,1]^ and is defined by Q. After that we will show that the closure of the 
outer bound can be achieved. 

An interesting observation is that the outer bound in Proposition |III.2| does not depend on 5i, 
i = S,R for qs < 5s and qR < 5r. Furthermore, if g^ is increased over 5i for i = S,R has no 
effect because the value of min(5j, g^) is bounded by 5i. 


In order to obtain the closure we have to solve two optimization problems. By replacing Xr 
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by X and A/j by y, the optimization problems are 

q's(l — qs) [PsD + (1 — Psd)Psr] Prd 


[Pi] max x = 

QS 


qs [PsD + (1 — Psd)Psr] 


and 


(1 — qs)PRD + ^' 5(1 — Psd)PsR (1 — qs)PRD + — Psd)PsR 


subject to y + “ qs)PRD 

{qs,qR) e [0,5s] X [o,5r] 

{x,y) e [ 0 ,1]2 


[P 2 ] max y = qnPnD 

IR 


(1 — Psd)Psr 


tX 


qRPRD 


(30) 

(31) 

(32) 

(33) 


[PsD + (1 — Psd)Psr] (1 — qR) [PsD + (1 — Psd)Psr] 


X (34) 


subject to X < qsil - qR) [Psd + (1 - Psd)Psr] 
{qs,qR) e [0,5s] X [o,5r] 

{x,y) e [0,1]2 

In order to solve [P 2 ], the differentiation of y with respect to gives 

dy 


dqn 


= P 


RD 


X 


(1 ~ qR)‘^ [Psd + (1 — Psd)Psr] 


The second derivative is negative since 

d^y 2 Prdx 


dq 


R 


(1 ~ qR.)^ [Psd + (1 — Psd)Psr] 


< 0 , 


(35) 

(36) 

(37) 


(38) 


(39) 


thus, the objective function is concave with respect to y. Solving the equation ^ = 0, we 
obtain the maximizing q\ where 


4 = 1 


X 


Psd + (1 — Psd)Psr' 


(40) 



17 


the corresponding maximum value of the objective function is 


y* = Prd - 2Pj 


RD 


X 


+ 


X 


PsD + (1 — Psd)Psr Psd + (1 — Psd)Psr 
Suppose that G (0,5j?), then 


[Prd — (1 — Psd)Psr\ ■ 

(41) 


(1 ~ ^rY [Psd + (1 ~ Psd)Psr] < x < [Psd + (1 ~ Psd)Psr\ , 


but since x <qs{l-qR) [Psd + (1 - Psd)Psr\ and qs e [0, 6s] then 


(42) 


X < 6l [Psd + (1 — Psd)Psr] ■ 


(43) 


We need to find the intersection of (42) and (43). If < 1 the intersection is an empty 

set, if 55 + > 1 then 


(1 — 6 rY [Psd + (1 — Psd)Psr] < x < 6^ [Psd + (1 — Psd)Psr] ■ 


(44) 


Suppose that the = 0 or = 5r, which is the case that x lies outside (44). If x lies 


outside (1441) on the right side then ^ is always non-positive and y is a. non-increasing function 


dqn 


of qR, thus = 0 and y* < 0. 

On the other hand if x is on the left side x < (1 — 6rY [Psd + (1 — Psd)Psr], then ^ 
is always non-negative and thus y is a non-decreasing function qR, hence q*ji = 6r and the 
maximum objective function is 

(1 — Psd)Psr 


y* = 6rPi 


-X 


6rPrd 


tX, 


(45) 


Psd + (1 — Psd)Psr (1 — 6r) [Psd + (1 — Psd)Psr\ 
for a; < (1 — 5rY [Psd + (1 — Psd)Psr] ■ (46) 


The initial constraint of the [P 2 ] should also met for q}^ = 6r, we have an additional condition 

X < 55(1 — 6r) [Psd + (1 — Psd)Psr]- 
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Summarizing, the closure obtained by [P 2 ], If < 1 then 

(1 — Psd)Psr __ SrPrd 


y* = SrP, 


RD 


-X 


for 


PsD + (1 — Psd)Psr (1 — Sr) [Psd + (1 — Psd)Psr] 


X < ^ 5(1 — Sr) [Psd + (1 ~ Psd)Psr\ ■ 


X 


(47) 


(48) 


If ^5 + > 1 then 


P, 


y = 


RD 2Prd^ Pg^+(i_%^)Pg^ + PsD+i^-PsD)PsR Psd)Psr] if Xi < a; < a;2 

if X < Xi 


X p_ {^-Psd)Psr ^ _ 

SrPrd Psd+{i-Psd)Psr 


SrPr 


{^-Sr)[PsD + {^-Psd)Psr] 

where xi = (1 — Sr)'^ [Psd + (1 ~ Psd)Psr] and X 2 = (5| [P^p + (1 — Psd)Psr]- 

Following the same methodology we obtain the solution to the [Pi]. If ^5 + < 1 then the 

closure is 


X 


(1 — Psd)Psrx 


+ 


y 


Ss [Psd + (1 — Psd)Psr\ (1 — Ss)Prd [Psd + (1 — Psd)Psr] (1 — Ss)Prd 

PsD+a-Ps^D)%R^ + ^ ~ 


= 1, (49) 


If Ss + Sr>1 then if Prp(1 - + y < PrdSr the closure is 


X 


Psr{1 - Psd)x ^ ^ ^ 


Psd + (1 — Psd)Psr V [Psd + (1 — Psd)Psr] Prd 

Psd+{i-%^d)Psr ^ + y- ~ closure is 

(1 - Psd)Psrx , y 


(50) 


X 


+ 


= 1. (51) 


Ss [Psd + (1 — Psd)Psr] (1 — Ss)Prd [Psd + (1 — Psd)Psr] (1 — Ss)Prd 

In the previous sections, we obtained the closure of the outer bound of the stability region, 
which we show below that this closure is achievable. The sufficient conditions for stability for 


fixed transmission probabilities are given in Proposition III.l For qs ^ Ss and qr < Sr we have 


that the saturated throughput for the source and the relay are given by 


ys — 7s(l ~ Qr) [psD + (1 — Psd)psr] , 


(52) 
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= 9i?(l - qs)PRD- 


(53) 


From (52) we obtain that 


Qs = 


Ps 


— Qr) \psd + (1 — Psd)psr\ 


By substituting (54) into (53) we have 


Pr = Qr 


1 - 


Ps 


Prd- 


(54) 


(55) 


(1 ~ Qr) [psd + (1 — Psd)psr]. 

After replaeing /i|j with Xr, total = Xr + pso+(iL°Jsd)psr ^^ identieal 


with the expression in (34) whieh eorresponds to the outer bound. As a result, a point of the 
saturated throughput can be controlled to any point on the boundary of TZ which is given in 
Proposition III.2[ The same argument holds for fXg and the expression in (30). 

The previous concludes that the closure of the outer bound of the stability region is achievable. 


VI. Conclusions 

In this paper, we studied the effect of energy constraints on a relay-aided wireless network in 
which both source and relay have energy harvesting capabilities. The source and the relay nodes 
also have external arrivals and network-level cooperation is employed, i.e. the relay forwards a 
fraction of the source’s traffic to the destination. 

We derived necessary and sufficient conditions for stability of the above cooperative communi¬ 
cation scenario and we also obtained the exact maximum stable throughput region. Interestingly, 
the closure of the inner and the outer bound is identical. A key insight of this work with an impact 
on the design of relay-assisted networks with energy limitations is as follows: when the aggregate 
charging rate is above one and both the source and the relay lie in the intermediate traffic regime, 
then the system has identical performance with the network without energy constraints. Otherwise 
stated, in the above setting, the energy limitations are transparent to the network operation, which 
is also demonstrated by the non-linear behavior of the bound of the maximum stable throughput 


region. 
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This work provides a step in eonneeting information theory and networking by studying the 
stable throughput region metric. Additionally, it sheds light on the relatively unexplored and 
important domain of energy harvesting and assesses the effect of that on this important measure. 

Future work will include the characterization of the stable throughput region using multi-packet 
reception instead of the erasure channel with collisions. 
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